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Abstract. Stimulated Brillouin scattering in optical waveguides is a fundamental 
interaction between light and acoustic waves mediated by electrostriction and 
photoelasticity. In this paper, we revisit the usual theory of this inelastic scattering 
process to get a joint system in which the acoustic wave is strongly coupled to the 
interference pattern between the optical waves. We show in particular that, when 
the optoacoustic coupling rate is comparable to the phonon damping rate, the system 
enters in the strong coupling regime, giving rise to avoided crossing of the dispersion 
curve and Rabi-like splitting. We further find that optoacoustic Rabi splitting could in 
principle be observed using backward stimulated Brillouin scattering in sub-wavelength 
diameter tapered optical fibers with moderate peak pump power. 


1. Introduction 

In the last few years, the held of cavity and Brillouin optomechanics, a branch of physics 
which focusses on the interaction between light and tiny mechanical and acoustical 
resonators, has drawn widespread interest because of key fundamental observations such 
as resolved-sideband cooling, optomechanically induced transparency, quantum coherent 
coupling, or Rabi-like oscillations PEI Eli]. Rabi oscillations were initially described 
as damped periodic oscillations of an excited atom coupled to an electromagnetic cavity 
in which the atom alternately emits and reabsorbs photons [5] . From a theoretical point 
of view, this remarkable and ubiquitous phenomenon can be readily described as a joint 
system of two-coupled oscillators. If both oscillators are uncoupled they both share the 
same degenerate eigenfrequency. But when they are strongly coupled, the degeneracy is 
removed and the frequency is split in two distinguishable eigenfrequencies corresponding 
to the odd and even supermodes of the joint system. The frequency difference between 
the two eigenfrequencies is called the Rabi frequency m- When the same effect occurs 
at the quantum level, meaning that the two oscillators are for instance an atom and 
a photon embedded in a cavity, the new eigensolutions are coherent superpositions of 
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the two particles. When the two particles are not interacting with each other, they 
have their own dispersion relations ki{uj) that may cross each other. However, if the 
strong coupling regime is achieved, the joint system dispersion curve displays an avoided 
crossing and Rabi splitting. Since its hrst observation, this concept has been extended 
to cavity optomechanics where mechanical and optical modes of high quality-factor 
resonators can be strongly coupled, and to many other joint systems including photons, 
phonons, excitons, polaritons or plasmons |2lElEl0i. 

In this work, we show that strong coupling regime and its formalism can also be 
applied to another strong photon-phonon interaction, which is known as stimulated 
Brillouin scattering (SBS). Specihcally, we theoretically demonstrate that such a strong 
coupling regime can be achieved between the acoustic phonon and the optical pump and 
Stokes waves using backward SBS in a subwavelength-diameter tapered optical hber 
P Hn]. Associated with the strong coupling is the observation of an avoided crossing 
and Rabi-like splitting between the optical and acoustical dispersion branches. 

The paper is organized as follows: First, we describe the principle and methodology 
of the stimulated photon-phonon interaction. Then we derive from the standard theory 
of SBS a joint two-level system, as commonly used in the Rabi problem. We investigate 
from this system the specihc conditions that allow for the achievement of strong 
optoacoustic coupling regime in tapered optical hbers. Both silica and chalcogenide 
glass materials are compared with different taper diameters and optical pump powers. 
Finally, we demonstrate that the strong coupling regime and Rabi-like splitting could be 
observed in a sub-micron (0.7pm) chalcogenide hber taper with moderate peak pump 
power (80 W). 

2. Principle and Methodology 

Let us hrst describe the stimulated photon-phonon interaction under consideration. 
When coherent laser light is coupled and guided into an ultrathin optical hber, as that 
shown schematically in Fig. [T| light generates and interacts with several types of acoustic 
waves mm- Here we assume only one acoustic pressure wave and we will hnd below 
the specihc conditions to do so in a hber taper. As scketched in Fig. [TJ SBS is an inelastic 
scattering whereby two frequency-detuned optical pump (in red) and Stokes (in blue) 
waves coherently interact in a dielectric material, giving rise to an optical interference 
pattern (purple) which generates an acoustic wave (green) from the electrostrictive force. 
Simultaneously, the photoelastic ehect creates an index grating that travels at the speed 
of hypersound (a few thousand of m.s“^ in silica and chalcogenide glasses). Pump light 
then scatters in the backward direction on the index grating that acts as a Bragg mirror, 
with a Doppler downshift that corresponds to the frequency detuning between the two 
optical waves. When the frequency detuning is equal to the phonon frequency, phase¬ 
matching is achieved and thus one gets amplihcation of the Stokes signal. In the weak 
coupling regime, both the optical interference pattern contrast and acoustic wave grow 
along the optical waveguide, whereas in the strong coupling regime, they hybridized 
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Figure 1. Principle of strong coupling regime using backward stimulated Brillouin 
scattering in an optical fiber taper. A strong pump wave (red) coherently interacts 
with a counterpropagating probe wave, giving rise to an optical interference pattern 
which generates and couples to an acoustic wave. Under strong coupling regime, the 
interference pattern and the phonon hybridized and alternately exchange energy over 
period of few microns. 


and alternately exchange energy over short period of a few microns, associated with the 
phonon lifetime. 

From a qnantnm point of view, a pnmp photon is annihilated to create both a 
Stokes photon and an aconstic phonon, as depicted in bottom left of Fig. [TJ Similarly 
to the Rabi problem, the two pnmp and Stokes photons can be seen as the two levels 
of an artihcial atom. In onr case however, the two conpled objects are actnally the 
phonon and the optical beating between the pnmp and Stokes waves which are conpled 
by the electrostrictive force. As a result, our joint system is fundamentally different 
from the exciton polariton picture. In the exciton polariton case, the two level system 
is supposed to have only one particle. It is then described as an exciton with a given 
resonant frequency, as sketched by the spring in bottom left of Fig. The photon 
interacting with the exciton is in fact a signal probe tuned at the resonant frequency. The 
experiment usually consists in tuning the probe signal by adjusting its angle (wavevector) 
or changing the temperature to affect the resonant frequency of the exciton. In the SBS 
case, the pump and Stokes light replace the two levels since up > us (see bottom right of 
Fig.ig. However, those two levels are not bounded by a resonant frequency as dehned by 
the gap of a two-level atom. In SBS, the Brillouin frequency shift (BFS) is given by the 
phase-matching condition between the acoustic properties of the mechanical structure 
formed by the hber taper and the optical waves. The resonant nature of the process 
mostly depends on the acoustic properties of the hber. This is illustrated in bottom 
right of Fig. where the phonon plays a role of spring. Since the pump and Stokes 
waves are coherent laser lights that are conveniently produced and manipulated, we 
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suggest to use them for probing the coupling with the acoustic phonon. 

To observe Rabi-like splitting in such a two-level system, it is also required to limit 
the optoacoustic interaction to ideally one guided acoustic mode, which is not possible 
using a bulky standard optical hber with a 125-/im cladding. This can however be 
achieved using an optical hber tapered down to a sub-micrometer diameter hber, to 
get only a few discrete acoustic modes. To that end, we hrst solved the dispersion 
equation of acoustic modes /3a (^) in a silica rod as a function of the taper diameter. 
This dispersion relation reads as mi 

— (9^ + /^DJi (H Ji (9a) - (H Ji (9a) = 4 / 3 aPgJi [pa] Jo (#)) 

a 

with p = ypl — (3^ and q = — jd'l where Vl and Vt are the longitudinal (P) and 

shear (S) acoustic velocities, respectively. Ji are the Bessel functions, a is the hber 
taper radius. We then did the same work for the optical modes, as was done in Ref. [12] . 
The optical dispersion equation for a step index hber of core refractive index ni and 
cladding refractive index n2 is 

[ Jm( 7 ia) nl K'^{i 2 a) \ m / 1 f ^ f 2 ) 

\ 7 iaJm( 7 ia) ’^2720X01(720)] \7i 72 /V 7 i nl'ylj' 

where 71 = 72 = \j(3“^ — and denotes the modihed function of 

the second kind with the prime denoting diherentiation with respect to the argument. 
ko is the propagating wave vector in vacuum. For each integer value m, the eigenvalue 
/3 is the ehective propagation constant along the hber axis of the given mode. This 
leads to the propagating constants of the pump (P) and Stokes (S) waves /3p(a;), and 
Then the phase-matching condition (3p{u + AfJ) — (dsi^) = /3a(Aff) sets the 
detuning frequency AfJ at which Brillouin scattering occurs. Figure [^typically shows a 
numerical simulation of the acoustic wave spectrum generated in a silica hber taper, as 
a function of frequency detuning Af2 and of the hber diameter [TTl ITU] . The color plot 
that represents the overlap between acoustic and optical modes actually corresponds 
to the scattering efficiency. As can be seen, there are several types of acoustic waves 
including longitudinal, shear, and surface waves around 5 GHz ra. We can also see 
several avoided crossings in Fig. that appear when a pressure (P) wave (dotted) 
crosses a shear (S) wave (dashed). Their polarizations are actually no longer orthogonal 
and thus they can interact with each other, giving rise to such avoided crossing. We 
must stress here that they are not related with the strong coupling regime described in 
the next section. Note also that around 5 GHz there is no anticrossing between the two 
surface Rayleigh waves because they have orthogonal polarizations. The squared red 
area in Fig. shows that, for sub-wavelength diameter hber, there is a small frequency 
range where we can isolate a single acoustic mode around 8 GHz. We will investigate 
the strong coupling regime in this area, which is sufficiently far from avoided crossing. 

Additionally, to reach the strong coupling regime, it is also needed to have a 
sufficiently strong coupling strength to transfer the energy back and forth between the 
optical beating and the phonon. As a result, large optical pump power is required. 
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We will see thereafter that chalcogenide glass-based fiber tapers with large Brillouin 
gain are the best potential candidates enabling the strong conpling regime. Also, it is 
important to stress here that we will not consider the radiation pressnre in onr model 
becanse it is negligible for the optical fiber tapers nnder investigation. Recent works 
have however shown that Brillonin gain conld be significantly enhanced in nanoscale 
and large step index wavegnides when radiation pressnre enters into play or becomes 
comparable to the electrostrictive force PSIIII]. In onr model, we will also neglect the 
anti-Stokes scattering and the second-order Stokes scattering [15]. That assnmption will 
be discnssed later in [Appendix E 
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Figure 2. Color plot of the coupling efficiency between the acoustic modes with the 
fundamental optical mode of a silica fiber taper as a function of the acoustic frequency 
(horizontal axis) and of the taper diameter (vertical axis). Several avoided crossings 
appear when a pressure (P) wave (dotted) crosses a shear (S) wave (dashed). The 
squared red area isolates a single acoustic mode, required for the strong coupling 
regime. 



3. Theory 

In the following, we reformnlate the nsnal theory of SBS, initially derived by Boyd [T5] . 
to yield a joint system that can experience Rabi splitting. The pnmp (P) and Stokes 
{S) fields and the material density p can be expressed as, 

Pp(x, t/, t) = Ap{z, t)Ep{x, + C.C., 

Es{x, y, z, t) = As{z, t)Es{x, y)e<f^sz-i^st) ^ ^ 

p{x, y, z, t) = po + [p{z, t)n{x, y)R(daz-nat) | _ 


( 3 ) 

( 4 ) 

( 5 ) 
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where po is the mean density of the medium, (3p^s,a and ojp^s and are the wavevectors 
and frequencies and acoustic wave, respectively. Ap^g ^ are slow varying 

enveloppes and Ep^g and TZ are the transverse profiles of each wave. Note that the P- 
wave model for an acoustic wave in a submicron fiber taper is appropriate if the frequency 
detuning range surrounds a P-wave branch and is far enough from the anticrossing with 
S-wave branch, as shown in the squared area in Fig. In this area, we can readily 
assume that only one acoustic wave can be possibly coupled to the two optical waves. 
Introducing these equations into the wave equation and omitting the complex conjugate, 
we find [I5] 


dAp djddAp lojp'ye 1 ~t 
dz doj dt 2?7,pcpo-Defr^ 

dAg ^ iujg'ye 1 

dz du dt 2ngcpo Des^ 


( 6 ) 

(7) 


where ^ jg electrostrictive constant, np, ng are the refractive indexes for 

the Pump and Stokes waves respectively, c is the speed of light in vacuum, Dgff is a 
factor coming from normalization that is discussed in Appendix A Note that this set of 
equations describes the nonlinear propagation of P and S waves but it does not exhibit 
the optical interference pattern. To that end, we need to combine them in order to get 
the equation of the optical beating. 


3.1. Poynting vector: Optical beating 

Since the optical beating is a local intensity modulation induced by the superposition 
of P and S waves, it can be described by the z component of the Poynting vector 

S = E xH, (8) 

= Sp + Ss + Ep X Hg + Es X Hp, (9) 

—¥ 

where S p,s are the Poynting vectors for the Pump and Stokes wave. The two first terms 
of the right-hand side (RHS) describe the independent S and P-wave Poynting vectors. 
They are linked to power evolution by the flux $. The two last terms in the RHS 
displays the coherent local intensity modulations due to the interference. Those terms 
do not describe a power flow since they vanish after integration, however they describe 
a local intensity fluctuation, i.e., the optical beat note. Using the expressions of the 
Pump and Stokes waves, @ and Q, we can readily see that this intensity fluctuations 
varies as ^ where Af3 = jdp — {dg and Acu = ojp — ujg. The optical beating thus 

travels at a speed defined by = ^. It is around a few thousands m/s, the same speed 
of the acoustic wave, if phase-matching is satisfied. In the following we will transform 
equations ([^ and Q to derive an equation of the optical beating and its relationship 
with the phonon. Careful attention was paid to the omitted complex conjugate terms, 
however all the terms of equation (|^ are not detailed and the focus is put only on 
phased-matched terms. First, we introduce the amplitude of the magnetic field Bp. 
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In the plane wave approximation, that gives for the P-wave Bp{z,t) = —-^^^Ap[zR) 
Equations (|^ and Q yield 

(3p iup'ye 1 


dBp ^ dfd dBp 


dz 

dBs 


doj dt 

dp dBs 


Ps 2npcpo -Deff 

Ps iws7e 


pBs, 


dz du dt 


Pp 2nscpo Des 


P*B, 


( 10 ) 

( 11 ) 


Multiplying equation by B*s and the conjugate of equation (11) hy Ap, we find 


dAt 

dz 


B\ 


n dB*s 
A, ^ 


~ _ 

^ dcj dt ^ 

_d_i~ d^_ 
du dt 


lUpye 


2npcpo Heff 

Ps 1 


pAsB*s, 

pApB*p. 


dz du ^ dt Pp 2nscpo Dgff' 

Since P and S lights are propagating in forward and backward directions and they are 
only shifted in frequency by a few GHz, we can make the assumptions ^ ~ — 1 and 
n = ns — Up. Taking the sum of above equations, and neglecting the group velocity as 
explained in Appendix B we readily find 

dApB*s 


lUpUs'^eP 


ApB*p 


AsB*s 


( 12 ) 

dz ncpoDes \ Up ) 

On the RHS of equation we see that the optical beat note is clearly coupled with the 
phonon p. The term within brackets is usually referred as the “population inversion” 
in a two-level system problem. Indeed, if both waves were propagating in the same 
direction, each term would be proportional to the power Pk of the related k- 

wave, and be denoted the “population” of level k. However in the specific case of 
counter propagating waves, the name “population inversion” is not appropriate. Since 
each “population” term is proportional to its Poynting vector, its sign depends on the 
direction of the power flow along the z axis. If the S-wave propagates in the backward 
direction (SBS), the flux of the Poynting vector for the Stokes wave is negative leading 
to — _Es_ As a result, the term in brackets of the RHS of the equation should be 

named “total population”. This situation is very specific to the study presented here, 
and that it does not occur in co-propagating or static cases like the well known atom- 
photon coupling. In the atom-photon picture, the conservation law stipulates that the 
sum of the populations is constant, pu + P 22 = 0[l5]. However, for counterpropagating 
waves P and S, the Poynting picture shows that the conservation turns to <hp — $5 = 0 
that is unusual. 

In order to build a joint-system where the optical beating is coupled with the 
phonon, it is important that the coupling rate does not vary, in other words, to conserve a 
constant “total population”. To fulfill this criterion, we have to assume that anti-Stokes 
and second-order Stokes scattering remain negligible (see Appendix E). Under this 


assumption, we can multiply equation ( 12 ) with the previously omitted fast oscillating 
phase term add and subtract —iApEpHg{z,t), simplify with equationand 

























SBS revisited: Strong coupling regime and Rabi splitting 
take the Fourier transform, we get 


dEpH*s{z, n) 
dz 


iApEpH*s{z,n) 


iupu:sle 

4ncpo^eff 



— (13) 

LVs ) 


with EpHg{z,t) = ApBg{z,t)e^^^^^~^‘^^\ and the hat denoting the Fourier transform. 
This equation actually describes the optical beat note in the Fourier domain, with 
being the frequency detuning between Pump and Stokes waves. If the RHS is neglected, 
meaning that there is no nonlinear interaction, this equation becomes the dispersion 
equation of the optical beating. We can then see that the optical beating travels with 
the wavevector Afd that is related to the speed v = 


3.2. Acoustic eguation 

To build the other half of the joint system, a similar equation must be derived for the 
acoustic phonon. Using a standard description of a pressure wave for the acoustic model, 
we get 


- 2ifl 


dp 

dt 


+ {nl inTp) p - 21 /?, 


,dp 


^"'dz 


lePl 1 


nc Heff 


ApB 


pJDg. 


(14) 


The hrst and third terms of the left-hand side (LHS) of the equation show that the 
phonon travels at the speed = %. The second term of the LHS is the signature of 
the mechanical stresses induced by the phonon. It is zero at the resonant frequency fla- 
The RHS of Eq. (14) provides the coupling rate with the optical beating. To write this 
equation, we have to suppose that the phase-matching for backward SBS is satished 
and therefore that the acoustic wavevector is twice the optical one, as [3a = ‘2(3p. Now, 
we should describe the whole dispersion curve and not restrict our study to the phase¬ 
matching case only. We will thus not make this assumption and let Papl) varying in 
the surrounding of frequency resonance (H — -C H^). The — — iUT^ term then 

reduces to —iUFs, which stands for the acoustic phonon losses. 

In order to see the dispersion equation of the phonon, we will rewrite the above 
equation by keeping the fast oscillating phase term ^ then add and subtract 

and straightforwardly get 


dPap + 




dt 


dz 


~ m{Q)p^ + ikpoEpHg, 


(15) 


where p^ 


p{z,t)&^^°-^ -I-c.c. and 


m(H) 


iUFs 

1 lePl 1 

2/3aVl nsc T>eff ’ 


(16) 

(17) 


represent the acoustic losses and the nonlinear coupling rate, respectively. The 
dispersion equation can be derived from the Fourier transform of equation Making 
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the assumption that the detuning from the Brillouin frequency shift is low -C fla, 

we hnd 


~ = m{Vt)p{z, VL) + ikpoEpH^{z, Q), (18) 

where /3a(^) = /3a + (^ — is the dispersion relation of the phonon and p(z, Q) is 

the Fourier transform of p^. 


3.3. The joint system 


as a joint system with 2EpHg denoting the optical 

beating 


We can write equations 


13 


and 


18 


/ 2EpHl \ _ . / koo kop 
dz \ p j Y kpQ kpp 



(19) 


where koo = ^/3 and kpp = (3a{El) + u{fl) are the eigenvalues of the uncoupled system. 
Those are the wavevectors of the optical beating and of the phonon, respectively. 
Note that in the expression of kpp, we have u{Ua) = that is related to the 
Brillouin linewidth F^, and thus the phonon decay rate or lifetime. Since this term is 
imaginary, it gives an imaginary part to kpp that stands for the phonon propagation 
losses. The coupling coefficient are kpo, previously dehned in equation 0 , and 
kop = — 2nopoD’'a eigenvalues K± of the joint system are solutions of 

characteristic polynomial 


koo) (-^i kpp) kopkpo 0 . 

If kopkpo = 0, the coupling rate is null and there are two independent eigenvalues with 
dispersion curves that cross each other. If the coupling term is non zero, it is more 
convenient to write the equation as the form 

(^koo T kpp )T kppkoo kopkpo 0, 


with 


k k — 

^op^po 


/doEpc go 


8nVa Aeff [uip UJs 


Pp ^ Ps 


( 22 ) 


where is fhe Brillouin gain dehned in Appendix C and A^ff = the effective 
area of the mode-held diameter (MFD). This equation possesses two solutions if the 
discriminant is positive 


^ (^oo T l^pp) d (^kppkoo kopkpo) 'R 0. (^ii) 

At the crossing of two dispersion curves (D = Dq), where the phase matching condition 
is fulhlled A/3 = (3p — /ds = (3a, neglecting the phonon losses leads to koo = kpp. We can 
rewrite this equation as 


A i^koo kpp) ~\~ dikopkpo ^ 0 . 


(24) 
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Here we see that, at the crossing point of the uncoupled dispersion curves, where the 
real part of koo equals the real part of kpp, only the coupling term and the imaginary 
parts of koo and kpp still remains in the equation. Therefore, if the sum of the remaining 
terms is positive the eigenvalues are 

K± = ]^ [koo + kpp ± a/a| . (25) 

We hrst compute the result of this equation, neglecting the phonon losses u{fl) = 0. 
This is illustrated in Fig.j^that shows the wavevector K± as a function of detuning for 
a 0.7 yum mode-held diameter chalcogenide tapered optical hber, with 100 W and 1 mW 
pump and Stokes power. The Rabi splitting is clearly visible in Fig. as an avoided 



Figure 3. Dispersion curves K±{^) showing acoustic Rabi splitting and avoided 
crossing in a 0.7 /rm mode-field diameter As 2 Se 3 tapered optical fiber with 100 W 
pump power and 1 mW Stokes power. The parameters are from table [Dll of [Appendix] 
l^and the phonon losses were neglected [20]. Dashed gray lines show the intersecting 
dispersion curves when there is no coupling. 

crossing between the optical beating and phonon branches. Specihcally, the splitting 
is characterized by a strong gap between the two curves of a/A. As a comparison, the 
dashed gray lines show the intersecting dispersion curves when there is no coupling. Note 
that for the simulation shown here, we considered an acoustic velocity of 2250 m.s“^ 
for chalogenide glass [S], so the splitting occurs around 8 GHz. Now, we will take into 
account the phonon losses and see how it will affect the Rabi splitting. Furthermore, 
we will dehne a criterion for the strong coupling regime and study how this criterion 
evolves with the material nonlinearity, the pump power and the hber taper diameter. 
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In the last paragraph, the imaginary part of koo and kop were neglected making each 
wave lossless. In the exciton-polariton picture, the weak and strong coupling regime 
are usually separated by comparing the coupling rate g and the exciton-photon decay 
rates 'jph, and ■jex- If S' 3> jeT, Iph the system is in the strong coupling-regime whereas if 
g <C 7 e 3 ;, Iph the system is in the weak coupling regime [16]. Our case using SBS is very 
similar. If we do not neglected the phonon losses in (24), we get the following inequality 

a/A = 2^kopkpo > |m(0 )|. (26) 

where the right-hand side y/A plays a role similar to the coupling rate g, with the 
difference that our eigenvalues are wavevectors, and a/A is thus homogeneous to 
instead of The LHS term u{Vt) dehned by equation (16) is related to the phonon 
decay rate, yea, in the exciton-polariton picture. To discriminate the weak coupling from 
the strong coupling regime, it is convenient to rewrite equation (26) with the splitting 
If p -C 1, the joint system is in the weak coupling regime whereas if 


ratio p = 


p 1, it is in the strong optoacoustic coupling regime. 

Figure compares these two distinct regimes for an 0.7 /im-diameter optical hber 
taper based on either chacolgenide or silica glass materials, whose parameters are 
listed in Dl^ in Appendix D The chosen glass materials are silica and chalcogenide 
As 2 Se 3 glasses because they are available as micro and nanowires [TUI ITT] . This 
comparison allows us to show the signihcant role of the Brillouin gain in the splitting 
ratio. Figure l^a) shows in particular the case where p < 1, for a pump power of 5W, 
in the gray area of Fig. |^b). In this case, the nonlinear process is not strong enough 
and no splitting can be observed. Figure |^c) clearly shows instead the strong gap \/A 
between the two branches for a pump power of 80 W. Rabi splitting occurs when the 
gap then exceeds a certain number of Brillouin linewidths dehned by |'u(f2)|. 

This number of linewidths is indeed the splitting ratio p that is plotted in Fig. |^b) 
as a function of the pump power. One can see that, for a the silica-based taper, the 
Brillouin gain is not strong enough to induce any splitting even at high pump power, 
whereas the chalcogenide hber taper enables the splitting from a pump power just above 
10 W. This is simply due to the the fact that Brillouin gain in As 2 Se 3 hber taper is 
almost 200 times that of silica (See Table Dl). The splitting ratio is further illustrated 
in Fig. igc ) for 80 W power. In this case p = 3, we thus expect a gap between the two 
branches of 3 linewidths. 

A further comparison between Figs. [^ and shows that the optical beating branch 
vanishes. This is due to the fact that optical beating experiences negligible loss (1 
dB.m“^), and thus its branch is very thin and beyond image resolution. Unlike the 
optical beating, acoustic phonon experiences much larger loss, leading to a much wider 
branch than the optical one. Consequently, we can deduce that the resonant nature of 
the process is dehnitely related to the phonon lifetime. 

Once we have seen how the splitting ratio p evolves with the optical power, it is 
legitimate to investigate the influence of the hber taper diameter. In our model however. 
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(a) (b) (c) 



Detuning (GHz) Punnp Power (W) Detuning (GHz) 

Figure 4. (a) Dispersion curves K±{n) for 5 W Pump power and 1 mW Stokes power 
in a 0.7 gm mode-field diameter tapered As 2 Se 3 fiber, (b) Wavevector splitting ratio p 
As 2 Se 3 and silica fiber taper versus the pump power, p measures the gap between two 
split curves in number of linewidth. The gray area delimits the weak coupling regime, 
(c) Dispersion curve K±{n) for 80 W Pump power and 1 mW Stokes power. Optical 
loss is set to 1 dB.m“^. 

it is more convenient to take into account the mode field diameter (MFD). Using this 
assumption, we compute the splitting ratio for different MFDs and 50 W pump power. 
The results are summarized in Fig. |^a) for chalcogenide and silica fiber tapers. Once 
again, the silica fiber taper does not provide enough coupling strength even at submicron 
scale. This is mainly due to the low effective index of silica fiber taper, as shown in 
Fig. I^b). For a diameter close to 0.5/im the effective index of the fundamental mode 
approaches 1, meaning that light is not highly confined into such a low step-index fiber 
taper. On the other hand. Fig. |^b) shows that the chalcogenide fiber taper has a higher 
effective index, enabling waveguiding with diameter down to 0.4 fim and a wide range of 
mode-field diameter to observe the splitting. This is depicted in Fig. |^a) where splitting 
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Figure 5, (a) Splitting ratio p for different mode-field diameters, 50 W pump power 
and 1 mW Stokes power. As 2 Se 3 and silica fiber taper are compared. Gray area 
represents the weak coupling region, (b) Effective index calculation for a fiber taper 
in silica and As 2 Se 3 with respect to the taper diameter. 
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can be observed with 50 W pump power even with a 1 pm diameter. Moreover, such 
peak pump power are easily achievable using nanosecond pulses at 1.55 pm and Erbium- 
doped fiber amplifier (EDFA). As a consequence, we can expect that chalcogenide fiber 
tapers are best candidate for an experimental demonstration of strong SBS coupling 
regime. 

5. Conclusion 

In conclusion, we have revisited the theory of stimulated Brillouin scattering in optical 
waveguides and demonstrated that strong photon-phonon coupling regime can be 
achieved between the acoustic wave and the two optical waves, giving rise to avoided 
crossing and Rabi-like splitting. This has been demonstrated by deriving from the 
coupled equations of SBS a joint system that combines both the phonon and the optical 
interference pattern that results from the coherent superposition of the pump and Stokes 
waves. It has been shown that, when the nonlinear coupling rate becomes comparable 
to the phonon decay rate, the joint system enters in the strong coupling regime and 
exhibits two splitted eigenvectors K±{fl). This regime was further investigated in detail 
as a function of the Brillouin gain, the pump power, and the mode-field diameter. As 
a conclusion, our results show that the strong coupling regime could be observed using 
backward SBS in a chalcogenide-glass 0.7/im-diameter optical fiber taper, simply by 
tuning the probe frequency around the Brillouin frequency shift. A splitting of the 
Brillouin spectrum should be observed at high pump power. 
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Appendix A. Normalization 


The fields are normalized so that the modulus \2AB* 
Thus the effective length Des can be written as 


in Eq. (12) is the optical power. 


JJ\Es\^dxdy^JJ\R\Mxdy 
11'l^(x,y)EsE*pdxdy 



where a = is the waist and MFD is the mode field diameter (after Ref. |18]). We 
can also define the effective mode area as Ag// = 
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If we do not neglect the group velocity in Eq. (12), we get 

d /3 r dAp 


dApBl 
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dz dui \ dt 


-B*s - Ap^^ 
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ApB*p ApB 


s^s 
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Note that if the two P and S waves were propagating along the same direction, the left- 

giving rise to a forward optical beating 


hand side (LHS) would become 


dz 

travelling at the mean group velocity. In the counterpropagating case, it is much more 
difficult to dehne a group velocity for the optical interference pattern. However, we can 
show that the effective group velocity is somehow bounded and has a negligible effect 
on our main results. We can see in Fig. [^a) three different scenario. First, if the pump 
is a continuous wave, = 0, and the Stokes wave is an optical pulse, the intensity 
fluctuations due to the superposition of both waves is located at the same location as 
the Stokes pulse. The optical beat pattern thus follows the Stokes pulse and propagates 
at the same group velocity. It translates into the equation by the fact that = 0. 
Note that the interaction with the phonon is not likely to disturb this equilibrium if the 
pump is sufficiently strong compared to the Stokes wave |Ap| S> Fig. |^b) shows 
the opposite scenario where the Stokes is a continuous wave and the pump a pulse. In 
this case, the optical beating follows the pump pulse and its group velocity is therefore 
the pump group velocity. The striking consequence is that this group velocity is the 
opposite of the previous case, as Pump and Stokes propagate in opposite direction. From 
the equation, we can address a third interesting example ii Ap = Bf, the group velocity 
of the optical beating is then null. This is depicted in Figj^c), where both pump and 
Stokes waves are optical pulses and that the barycenter of both pulse is static. This 
last case is however very unlikely to happen since an increasing of Bg over time due 
to the phonon is clearly related to a depletion of the Pump Ap with much impact on 
the group velocity. As a result, if the group velocity of the optical interference pattern 

where the 


_ai ,dp 

du) ’ ' dui 


may be uncertained, it is certainly bounded by the relation A ^ 
limits correspond to the two hrst cases. More intriguingly, thought these boundaries 
are wide appart, they do not affect our main results. Figure B1 (d) shows a similar 
numerical simulation of the avoided crossing as in Fig. but with A = —^ in plain 
green and blue and ~ = +f5 with markers. There is no noticeable difference because 
T ~ ^ T is very small compared to A and thus negligible whatever the sign is. 
The second point we want to address is that in most Brillouin experiments, the group 
velocity is always positive since the pump power is greater than the Stokes one before 


depletion. Such assumption leads to 
velocity. 






dt 


leading to a positive group 
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(d) 



Figure Bl. (a) Pump in gray is a continuous wave and travels forward the z-axis, 
Stokes pulse light travels backward, (b) Pump pulse in black travels forward the z-axis, 
Stokes light in gray is continuous and travels backward, (c) Both lights are pulses and 
travel in opposite direction. The resulting wave packet is static, (d) Same computation 
as in Fig. 3 but with ^ in plain green and blue and with ^ + in markers 


Appendix C. Brillouin gain 


The Brillouin gain was defined using the method described in [19]. From the above 
equations ([^ and (0 , we find 
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These equations first show that the optical powers propagate at the same group velocity. 

Moreover, if we neglect the nonlinear coupling term and then substract both equations, 

—¥ 

we get the Poynting theorem for plane waves, div S= 0 and the conservation law 
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Using the phonon equation (14) in steady-state regime, in order to substitute p, we 
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and assuming cap ~ we thus obtain the usual power equation evolution of SBS, as in 
Ref. 119] 
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where 
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Appendix D. Glass parameters 

In onr simnlations, we nsed and compared silica-based to chalcogenide As 2 Se 3 glass- 
based hber tapers, becanse they are readily available and mannfactnred DU. Moreover, 
recent experiments have shown great potential of these glass materials for SBS 
applications PHO]. The table below snmmarizes the optical and acoustical parameters 
of these two materials. Note that the As 2 Se 3 Brillouin linewidth is taken larger than 
the one in [20] because of the spectral broadening due to hber coating [2]. 


Table Dl. Table of optical and acoustical parameters of silica-based and 
chalcogenide-based glass materials. 



silica“ 

As2Se3“ 

refractive index at 1.5 gm 

1.45 

2.8 

Brillouin frequency shift (GHz) 

11 

8 

Brillouin linewidth (MHz) 

30 

50 

Acoustic velocity (ms“^) 

5600 

2250 

Brillouin gain (mW“^) 

3.3.10-11 

600.10-11 

Mean mass density (Kg.m“^) 

2200 

4640 


“See [IS]; '’See |20|. 
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One of the main hypothesis for the strong coupling Brillouin regime was to consider 
both the anti-Stokes and second-order Stokes scattering as negligible. This assumption 
is valid in the weak pump to Stokes conversion regime. However, adding the anti-Stokes 
wave as a higher energy level in the system equation does not signihcantly change our 
results. As shown in Fig. El, anti-Stokes scattering involves an oncoming acoustic 
phonon A, that must not be mistaken for the forward propagating phonon B, involved 
in the joint pump-Stokes svstem[T5]. The scheme illustrates that anti-Stokes scattering 
can then be considered as an additional loss for the pump level. As a result, the model 
could be considered as a 3-level system [5]. The pump depletion due to this phenomenon 
would then be associated to a decay rate related to the room temperature phonon 
population. However, in optical hbers, anti-Stokes scattering is a much weaker than 
Stokes scattering due to the low phonon population and the fast phonon decay rate. 
Therefore, it is usually neglected and particularly in the stimulated regime [21]. In the 
strong coupling regime, the coupling must be stronger than the phonon decay that is 
around few microns. Over such a short distance, the pump decay rate associated with 
anti-Stokes scattering is thus negligible. Note that second-order Stokes scattering can 
be neglected for the same reasons. 


Anti-Stokes 



h 

Y*nmnmr- 
j Phonon A 

Pump 

Stokes 

J Phonon S 


Joint-System 


Figure El. Scheme of the anti-Stokes scattering process above the joint pump-Stokes 
system. 
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